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tower located at KSC is determined as a function of wind direction. 

The roughness length estimates, which were calculated with wind profile 
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determined for thirty-nine wind and temperature profiles. Most of the 
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TECHNICAL MEMORANDUM X-53690 


AN ANALYSIS OF THE ROUGHNESS LENGTH ASSOCIATED WITH 
THE NASA 150 -METER METEOROLOGICAL TOWER 

SUMMARY 


The surface roughness length at the NASA 150-meter meteorological 
tower located at KSC is determined as a function of wind direction. 

The roughness length estimates, which were calculated with wind profile 
laws consistent with the Monin and Obukhov similarity hypothesis, were 
determined for thirty-nine wind and temperature profiles. Most of the 
cases were obtained during the hours of 0700 and 1600 EST ? and the dura- 
tion of each test ranged between one-half to one hour. The mean wind 
speed data were obtained at the 18- and 30-meter levels, and the mean 
temperature data were obtained at the 18- and 60-meter levels. For 
those wind directions 0 in the ranges 0° ^ 0 < 150°, 180° - 0 < 240°, 
and 300° ^ 0 < 360° the roughness length is 0.23m; for those wind direc- 
tions in the ranges 150° i 0 < 180° and 240° ^ 0 < 300°, the roughness 
length has the values 0.51m and 0.65m, respectively. 


I. INTRODUCTION 


In recent years, low altitude winds and turbulence have been the 
object of extensive interest in the aerospace and meteorological com- 
munity for the design of space vehicles, buildings, bridges, antennae, 
aircraft, etc., turbulent diffusion, aircraft operations, and many 
other engineering and scientific problems. The National Aeronautics and 
Space Administration has addressed the problem of low level winds and 
turbulence in the context of space vehicle design. NASA personnel are 
now developing analytical models of launch vehicles which predict the 
response of these vehicles to various types of ground wind forcing func- 
tions. These forcing functions can be prescribed in terms of wind pro- 
files, discrete gusts, gust factors, and spectral estimates of turbulent 
wind fluctuations. 

To provide meaningful ground wind design data for these response 
and loading calculations, NASA has constructed a 150-meter meteorological 
tower at the Kennedy Space Center, Florida. This tower, described by 
Kaufman and Keene (1965), is located in the vicinity of Launch Complex 39 
and is situated in a well exposed area free of nearby structures which 
could interfere with the air flow. It is instrumented at the 18, 30, 



60, 90, 120, and 150 meter levels with Climet wind sensors (Model Cl-14) 
and Climet aspirated thermocouples (Model 016) located at the 18, 60, 
120, and 150 meter levels. 

To develop turbulence models based upon the data from this tower, 
the surface roughness length z Q of this site must be known. The param- 
eter is of no great interest in itself, but rather, its importance lies 
in the fact that it serves as a scaling length in the formulation of 
boundary layer wind profile laws based upon asymptotic similarity con- 
siderations (Blackadar, 1967) or heuristic considerations using eddy 
coefficients (Blackadar, 1965). In the absence of vertical velocity 
fluctuation data -- this is the situation at the NASA 150-meter meteoro- 
logical tower -- these profile laws permit us to calculate the surface 
fridtion velocity. The friction velocity in turn is used as a velocity 
scaling parameter which ultimately permits us to combine turbulence data 
in the form of spectra, cospectra, variances, etc., based upon similar- 
ity considerations. 


II. THEORETICAL BACKGROUND 


In the first 30 to 60 meters of the atmosphere, the Monin-Obukhov 
similarity hypothesis predicts that 


kz du 
u. dz 


0(z/L), 


( 1 ) 


where u is the mean 1 wind at height z, u* is the surface friction veloc- 
ity, k is von Karman 1 s constant with numerical value equal to approxi- 
mately 0.4, 0(z/L) is a universal function of z/L which is determined 
experimentally, and L is the Monin-Obukhov stability length given by 


L 


u^ c pT 

JL p 

kg H 


( 2 ) 


where Cp is the specific heat at constant pressure, p and T denote the 
mean density and Kelvin temperature, respectively, g is acceleration of 
gravity, and H is the vertical heat flux. The friction velocity u* is 
given by 

u* = J- u f w' , (3) 


i 

All means are time averages. 
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where u f and w 1 denote the longitudinal and vertical velocity fluctua- 
tions about the mean wind vector and the vertical heat flux is given by 


H = pc p w' T' , (4) 

where T f is the temperature fluctuation about the mean temperature and 
the overbar denotes the time averaging operator 

T 

( ) - i /< ) dt, 

o 


T being the interval of time over which the averaging process is per- 
formed. If 0(z/L) is known, it is possible to obtain the wind profile 
by integrating equation (1) and applying the boundary condition that u 
must vanish at z = z 0 , where z 0 is the surface roughness length. z Q 
may be interpreted as that height above the local mean level surface of 
the earth below which the flow has been disrupted to the extent that the 
flow is completely turbulent; i.e., the mean flow vanishes for z ^ z 0 . 

In terms of 0(z/L), the wind profile is given by 

u ' T { ln t ' I “‘I ' <5) 

-z 0 /L 

There are various ways to calculate the roughness length z Q from equa- 
tion (5); however, in any event, three pieces of information are required 
This information can be (1) the mean wind speed measured at three levels, 
or (2) the mean wind speed measured at two levels and either u* if H is 
unknown or H if u* is unknown, or (3) the mean wind speed measured at 
one level, u* and H. In the Monin-Obukhov theory u* and H are height 
invariant, so that a calculation of these quantities at more than one 
level gains little, unless the theory is incorrect. In that case, the 
wind profile law must be reformulated ab initio. However, the theory is 
most applicable near the surface of the earth so that the best estimate 
of u* and H is obtained at the lower levels near the surface of the earth 
say below 10 m. In the first case, one would evaluate equation (5) at 
three levels to yield three transcendental equations in three unknowns, 
namely, z Q , H and u*. These equations would be difficult to solve since 
it would require a tr ial-and-error approach on the computer in view of 
their complicated transcendental nature. In the case of the NASA 150- 
meter meteorological tower, we would have to use the mean wind speeds 
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obtained at the 18-, 30- , and 60-meter levels; however, it is question- 
able to employ the 60-meter level wind speeds since the Monin-Obukhov 
theory tends to fail at these levels. In the second case, we would 
evaluate equation (5) at two levels and solve for two unknowns, namely 
z Q and H if u* is known, and z Q and u^ if H is known. In the third case 
we would evaluate equation (5) at one level and solve for z Q directly. 
However, in both cases (2) and (3), we require w f data to calculate u* 
and/or H. At the present time, the NASA 150-meter meteorological tower 
does not possess the capability to obtain w 1 data; therefore, we must 
use technique (1) to calculate z Q or cast the theory in terms of param- 
eters that can be measured within the capabilities of the current instru 
mentation at KSC in order to avoid the problem of solving a set of 
complicated transcendental equations. The latter alternative is the 
most expedient approach and it will permit a calculation of z Q with wind 
data obtained at the 18- and 30-meter levels, levels at which the Monin- 
Obukhov theory is valid. 

The flux Richardson number is defined to be 


R 


f 



H 


PC P 


u 


2 d_U 

* dz 


( 6 ) 


Combining equations (1) and (6), we find 


^ = 0 (z/L) R f . 


( 7 ) 


The eddy heat conduction and viscosity coefficients are given by 


H 



d z 


( 8 ) 


and 


K 

m 


du 5 
dz 


( 9 ) 


where 0 is the mean potential temperature. 
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Upon combining equations (2) , (8), and (9), we find 


where 



L L 1 K 


m 


L' 


U * T 


du 

dz 



( 10 ) 


(U) 


Combining (6), (8) and (9), we obtain the additional relationship 


R 


f 



Ri, 


( 12 ) 


where Ri is the gradient Richardson number given by 


Ri = 


£. d Q 

T dz 
(— ) 2 ’ 


(13) 


Ri is called the gradient Richardson number since it is based upon the 
gradients of 9 and u. Most investigators postulate K^/l^ to be a func- 
tion of the gradient Richardson number. If this assumption is imposed, 
equations (7), (10) and (12) constitute three equations in four unknowns, 

namely, Rf, Ri, z/L and z/L' , and it is possible to express three of 
these unknowns in terms of the -fourth one. Thus, we can solve for Rf, 

Ri, and z/L as functions of z/L' and conclude from equation (1) that 




(14) 


where 0 x (z/L' ) is a universal function of z/L'. Integration of equation 
(14) yields a wind profile law identical in form to equation (5) with 
and L' replacing 0 and L, respectively, so that 
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u 


(15) 



t(z/L') 


where ^(z/L 1 ) is a universal function of z/L' given by 


-z/L' 

♦<z/L') - J — d| (16) 

-z 0 /L' 

provided 0 1 (z/L* ) is a known function. 

The derivative de/dz in equation (13) can be related to the vertical 
temperature gradient by differentiating Poisson’s law, 

0 = T( Pl /p) R/c P (17) 


and employing the condition for hydrostatic equilibrium. 


lE = 

dz 


" pg, 


(18) 


where Pi is the standard sea level pressure in appropriate units and R is 
the specific gas constant for air. The desired result is 


1 d0 = 1 /clT + 

G dz T \dz c J 

\ p/ 


(19) 


We have used the ideal gas law 

p = RpT (20) 

in deriving equation (19). Within the atmospheric boundary layer we can 
approximate 
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so that equation (13) reads 


Ri ^ 


£ 

T 


■dT 

_d 

du 

dz 


-£. 

2 



( 21 ) 


In practice, one usually obtains an estimate of Ri based upon the 
mean wind speed and temperature profiles and then calculates z/L 1 from 
an experimentally determined expression relating Ri and z/L 1 . Such an 
expression is predicted to exist according to this analysis. Ri as given 
by (21) is estimated with mean flow quantities that can be measured with 
data obtained from the NASA 150-meter meteorological tower. Upon deter- 
mining z/L 1 , one calculates \|/ with another experimentally determined 
relationship which is obtained by experimentally determining 0 1 (z/L l ) 
and then producing the integral (15) with the assumption that z Q /L f « z/L f 
so that the lower limit on this integral can be approximated as z 0 /L f ~ 0. 
The functions z/L' (Ri) and ^(z/L 1 ) have been determined by other investi- 
gators based upon simultaneous measurements of the wind profile, the 
longitudinal and vertical velocity fluctuations, and the mean temperature 
profile made at various tower sites around the globe. According to Lumley 
and Panofsky (1963), 


z 

L 1 


Ri 

(1 - 18R1) 1 / 4 


(Ri < -0.01) 


( 22 ) 


z 

L 1 


= Ri 


(-0.01 i Ri ^ 0.01) (23) 


z _ Ri 
L 1 1 - 7Ri 


(0.1 * Ri > 0.01). (24) 


Presumably for Ri >0.1 no simple relation exists between Ri and z/L 1 
since turbulence, if present at all, is relatively weak and the mean 
flow at various levels tends to be uncoupled. Those cases associated 
with Ri >0.1 were not examined in this study. The function \|/(z/L f ) 
that correspond to the expressions given by equations (23) and (24) 
are given by 



(-0.01 g Ri g 0.01) 


(25) 


i|f(z/L' ) = -4.5 

t(z/L') = ~7 fr (0.1 I Ri > 0.01). (26) 

Panofsky (1963) has graphically indicated the function ^(z/L* ) for 
Ri < -0.01, and this author finds that the function 

-7/1 1 

1,0674-. 0678 ln( n 

f /t » \ °- 01 

t(z/L') = .044 (Ri < -0.01) (27) 


faithfully reproduces his curve. This function is shown in figure 1. 

One should keep in mind that the relationships (22) through (27) have 
been deduced from measurements obtained at other sites. The Monin- 
Obukhov similarity hypothesis concerns the dependency of the gradients 
of 0 and u upon z/L and thus z/L 1 in view of our assumption concerning 
Kh/I^. This hypothesis is independent of the site, and it is only 
through the integration of equation (14) and application of the lower 
boundary condition (u = 0 at z = z Q ) that the site enters the picture. 

In short, kz u~ x du/dz does not explicitly depend upon the terrain 
features of the site, while u depends upon the site through z Q . Since 
the functional relationship between z/L 1 and Ri was derived directly 
from the Monin-Obukhov similarity hypothesis and since Ri is independent 
of the site, the function z/L* (Ri) is truly universal, so that it is 
reasonable to assume that equations (22) through (24) are also valid at 
KSC . Now, the function ^(z/L 1 ) is weakly dependent upon z Q through the 

lower limit on the integral representation of \j /, so that \j/ is site- 

dependent. However, if we replace this lower limit with zero the con- 
tribution to over the domain 0 ^ z/L ! l ^ 0 /L f is negligible, so that, 
although equations (25) through (27) are strictly valid for other sites, 
the dependence upon z G and thus the site is small. Thus, they can be 
safely used in studies involving KSC data. 

Upon determining the mean temperature and wind speeds at two levels 
it is possible to calculate Ri with equation (21). Based upon this esti 
mate of Ri, we can calculate z/L 1 with the aid of one of the equations 

(22) through (24). Finally, the function \|/(z/L f ) can be evaluated with 

the appropriate equation from the set (25) through (27), and thus u^ and 
z Q can be calculated with equation (15). 
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III. DATA PROCESSING PROCEDURES 


This analysis was based upon thirty-nine cases of turbulence at 
Kennedy Space Center. Most of these measurements were obtained during 
the hours of 0700 and 1600 EST, and the duration time of each test 
ranged between one-half to one hour. Mean wind speed data were obtained 
at the 18- and 30-meter levels, and mean temperature data were obtained 
at the 18- and 60-meter levels. Temperature at the 30-meter level was 
estimated by interpolating logarithmically between the 18- and 60-meter 
levels with the expression 


-'I In — 

T(z x )j- ff- , (28) 


where z ± , z^, and z 3 equal 18, 30, and 60 m, respectively. An esti- 
mate of the gradient Richardson number at the 23-meter level (geo- 
metric mean height between the 18- and 30-meter levels) was determined 
by assuming logarithmic distributions for the mean wind and temperature 
between these levels. The gradient Richardson number estimated in 
this manner is given by 


t(z 2 ) = t(z 3 ) - -mz 3 ) 


R1( V 


£ fT(zg) - T(z t ) + _£_1 

T V In <3&) V 

g z i 

r l 2 

U ( Z g) - u( Zl ) 

In (z 2 /z!) 


where 


z g ■ ^ 


(29) 


(30) 


and 


T 


2 In zp - (zp + 1) In z x 
z 2 - z ± 


■) f T(z 2 ) - T(2 X )1 
1 J \ In (zl/z!) J 


+ T(z x ) . 


( 31 ) 
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Then, Zg/L* was evaluated for each case by means of one of the three 
equations (22) through (24) corresponding to the appropriate Richardson 
number, L 1 was then assumed to be invariant with height, and ^(18/L 1 ) 
and ^(30/L 1 ) were estimated with one of the equations given by (25) 
through (27) corresponding to the appropriate Richardson number class. 
Equation (15) was then evaluated at the 18- and 30-meter levels to yield 
two equations in the two unknowns, u* and z Q , which are given by 

r u(z 2 ) - u( Zl ) *>. 

U* ■ k < V (32) 

\Ln(z 2 / z x ) - ^(z^/L 1 ) + \|f( z^/L*) J 


and 


z q = z 2 exp - j ^ ku ( z 2 ) + z 2 /L')^ , (33) 

where one first evaluates equation (32) for u* and then in turn evaluates 
equation (33) for z Q based upon this value of u*. 


IV. RESULTS 


The data used in the calculations in section III are shown in 
table I, and the surface roughness length as a function of wind direc- 
tion is shown in figure 2. To determine if there were any directional 
variations in zo> the data in figure 2 were averaged over 30° sectors 
beginning with 0° reckoned clockwise from north (see figure 3). The 
broken-line portion of this graph associated with wind directions 
between 240° and 270° was obtained by linearly interpolating between 
the results of sectors 210° through 240° and 270° through 300°. This 
diagram shows that in the sectors 150° through 180° and 240° through 
300° the roughness length is significantly higher than the roughness 
length associated with the other sectors. These high values of rough- 
ness length in these sectors can be attributed to the presence of trees 
upstream from the tower. The NASA 150-meter meteorological tower at KSC 
and the surrounding vegetation are discussed in a report by Kaufman and 
Keene (1965). 
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Based upon a sector average and the results shown in figure 3, 
table II shows the values of roughness length and the associated wind 
direction ranges appropriate for analyzing wind profiles from the NASA 
150-meter meteorological tower. The statistical error for the roughness 
length in sectors 1, 3, and 5 is .14m and is based upon 32 observations. 
Since the observations in the other sectors were not sufficiently 
numerous to obtain a reliable estimate of the statistical error, the 
result must be considered tentative. 


TABLE II 


Wind Direction (9) z Q 

Sector (degrees) (meters) 


1 

0° 

g 

0 

< 

150° 

0.23 

2 

150° 

g 

9 

< 

180° 

0.51 

3 

180° 

g 

9 

< 

240° 

0.23 

4 

240° 

g 

0 

< 

300° 

0.65 

5 

300° 

g 

0 

< 

360° 

0.23 
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FIG. 1. PANOFSKY’S + FOR Ri <-0.01 
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FIG. 2. CALCULATED ROUGHNESS LENGTHS VS WIND DIRECTION 
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